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The class of controlled synchronization systems under information constraints imposed by limited 
information capacity of the coupling channel is analyzed. It is shown that the framework proposed in 
A. L. Fradkov, B. Andrievsky, R. J. Evans, Physical Review E 73, 066209 (2006) is suitable not only 
for observer-based synchronization but also for controlled master-slave synchronization via commu- 
nication channel with limited information capacity. A simple first order coder-decoder scheme is 
proposed and a theoretical analysis for multi-dimensional master-slave systems represented in the 
Lurie form (linear part plus nonlinearity depending only on measurable outputs) is provided. An 
output feedback control law is proposed based on the Passification theorem. It is shown that the 
upper bound of the limit synchronization error is proportional to the upper bound of the trans- 
mission error. As a consequence, both upper and lower bounds of limit synchronization error are 
proportional to the maximum rate of the coupling signal and inversely proportional to the informa- 
tion transmission rate (channel capacity). The results are applied to controlled synchronization of 
two chaotic Chua systems coupled via a controller and a channel with limited capacity. 
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I. INTRODUCTION 

Synchronization of nonlinear systems, particularly 
chaotic systems has attracted the attention of many 
researchers for several decades [l|, [3]. During recent 
years interest in controlled synchronization has increased, 
partly driven by a growing interest in the application of 
control theory methods in physics [3, 0, |a, 0] . The de- 
sign of links interconnecting parts of cornplex systems to 
enable synchronization was studied in [3, 0, Q. These 
and related papers explore the possibility of modifica- 
tion of complex dynamical system behavior by means of 
feedback action. 

Modifying the behavior of complex interconnected sys- 
tems and networks has attracted considerable interest. 
The available results significantly depend on models of 
interconnection between nodes. In some works the in- 
terconnections are modeled as delay elements. However, 
the spatial separation between nodes means that model- 
ing connections via communication channels with limited 
capacity is more realistic. Hence analysis of the overall 
system should include both dynamical and information 
aspects. 

Recently the limitations of control under constraints 
imposed by a finite capacity information channel have 
been investigated in detail in the control theory litera- 
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ture, see [13, [HI, [13, [3, Q [S d and references and 
the references therein. It has been shown that stabi- 
lization under information constraints is possible if and 
only if the capacity of the information channel exceeds 
the entropy production of the system at the equilibrium 
[l3,[l3,[I3|- I'^ [iSjliMl ^ general statement was proposed, 
claiming that the difference between the entropies of the 
open loop and the closed loop systems cannot exceed the 
information introduced by the controller, including the 
transmission rate of the information channel. However, 
results of the previous works on control systems analy- 
sis under information constraints do not apply to syn- 
chronization systems since in a synchronization problem 
trajectories in the phase space converge to a set (a man- 
ifold) rather than to a point, i.e. the problem cannot be 
reduced to simple stabilization. 



The first results on synchronization under information 
constraints were presented in [19|, where the so called 
observer-based synchronization scheme 0, [20| was con- 
sidered. In this paper we extend this work and analyze a 
controlled synchronization scheme for two nonlinear sys- 
tems. A major difficulty with the controlled synchroniza- 
tion problem arises because the coupling is implemented 
in a restricted manner via the control signal. Key tools 
used to solve the problem are quadratic Lyapunov func- 
tions and passification methods borrowed from control 
theory. To reduce technicalities we restrict our analysis 
by Lurie systems (linear part plus nonlinearity depending 
only on measurable outputs). 



II. DESCRIPTION OF CONTROLLED 
SYNCHRONIZATION SCHEME 

Consider two identical dynamical systems modeled in 
Lurie form (i.e. the right hand sides are split into a linear 
part and a nonlinear part which depends only on the mea- 
surable outputs) . Let one of the systems be controlled by 
a scalar control function u{t) whose action is restricted by 
a vector of control efficiencies B. The controlled system 
model is as follows: 

x{t) = Ax{t) + B(p{yi), yi{t) = Cx{t), (1) 

i(i) = Az{t) + B^{y2) + Bu, y2{t) - Cz{t), (2) 

where x{t), z(t) are n-dimensional (column) vectors of 
state variables; yi{t), y2{t) are scalar output variables; 
A is an {n x n)-matrix; _B is n x 1 (column) matrix; C 
is an 1 X n (row) matrix, ip{y) is a continuous nonlin- 
earity, acting in the span of control; vectors i, i stand 
for time-derivatives of x(t), z(t) correspondingly. System 
([1]) is called master (leader) system, while the controlled 
system ^ is slave system (follower) . Our goal is to evalu- 
ate limitations imposed on the synchronization precision 
by limited the transmission rate between the systems. 
The intermediate problem is to find a control function 
U(-) depending on measurable variables such that the 
synchronization error e(t), where e(t) = z(t) — x(t) be- 
comes small as t becomes large. We are also interested 
in the value of output synchronization error e(t) = y2(t) 
~yi{t) = Ce(t). 

A key difficulty arises because the output of the master 
system is not available directly but only through a com- 
munication channel with limited capacity. This means 
that the signal yi(t) must be coded at the transmitter 
side and codewords are then transmitted with only a fi- 
nite number of symbols per second thus introducing some 
error. We assume that the observed signal yi(t) is coded 
with symbols from a finite alphabet at discrete sampling 
time instants tk = kTs, k — 0, 1, 2, ... , where Tg is the 
sampling time. Let the coded symbol yi[k] = yi(tk) be 
transmitted over a digital communication channel with a 
finite capacity. To simplify the analysis, we assume that 
the observations are not corrupted by observation noise; 
transmission delay and transmission channel distortions 
may be neglected. Therefore, the discrete communica- 
tion channel with sampling period Ts is considered, but 
it is assumed that the coded symbols are available at the 
receiver side at the same sampling instant tk = kTs, as 
they are generated by the coder. Assume that zero-order 
extrapolation is used to convert the digital sequence yi [k] 
to the continuous-time input of the response system j/i (t) , 
namely, that yi(t) = yi[k] as kTs < t < (k + l)^^- Then 
the transmission error is defined as follows: 



Sy(t)^yi(t)~yi(t). 



(3) 



On the receiver side the signal is decoded introducing 
additional error and the controller can use only the signal 
yi(t) = yi(t) + 6y(t) instead of y(t). A block diagram of 
the system is shown in Fig. [TJ 
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FIG. 1: Block-diagram for master-slave controlled synchro- 
nization (master system output j/i is transmitted over the 
channel) . 



We restrict consideration to simple control functions 
in the form of static linear feedback 



u(t) = ~Ke(t), 



(4) 



where e(t) — y2{t) —yi(t) denotes an output synchro- 
nization error, K is a scalar controller gain. The prob- 
lem of finding static output feedback even for linear 
systems is one of the classical problems of control the- 
ory. Although substantial effort has been devoted to 
its solution and various necessary and sufficient condi- 
tions for stabilizability by static output feedback have 
been obtained, most existing conditions are not testable 
practically [2l|, [22|. Since we are dealing with a non- 
linear problem further complicated by information con- 
straints, we restrict our attention to sufficient condi- 
tions for solvability of the problem and evaluate upper 
bounds for synchronization error. To this end we intro- 
duce an upper bound on the limit synchronization error 
where the supremum is taken over 

In the next two sec- 
tions the coding and decoding procedures are described 
and a bound on admissible transmission errors Sy(t) is 
evaluated. 



III. CODING PROCEDURES 

In the paper [l9| the properties of observer-based 
synchronization for Lurie systems over a limited data 
rate communication channel with one-step memory time- 
varying coder is studied. It is shown that an upper bound 
on the limit synchronization error is proportional to a cer- 
tain upper bound on the transmission error. Under the 
assumption that a sampling time may be properly chosen, 
optimality of the binary coding in the sense of demanded 



Q = sup lim ||e(i)| 

t — ^oo 

all admissible transmission errors 



transmission rate is established, and the relationship be- 
tween synchronization accuracy and an optimal sampling 
time is found. On the basis of these results, the present 
paper deals with a binary coding procedure. 

At first, we introduce the mcmoryless (static) binary 
coder to be a discretized map (7 : M ^ M as 



M[k] = (Mo - Mo,)/ + M, 



0,1, 



(8) 



q{y,M) = Msign{y), 



(5) 



where sign(-) is the signum function: sign(y) = 1, if y > 
0, sign(j/) = — 1, if y < 0; parameter M may be referred 
to as a coder range. Evidently, \y — q(y, M)\ < M for all 
y such that y : \y\ < 2M. Notice that for a binary coder 
each codeword symbol contains one bit of information. 
The discretized output of the considered coder is given 
asy = q{y, M) and we assume that the coder and decoder 
make decisions based on the same information. 

The static coder (l5|) is a subset of the class of time- 
varying coders with memory, see e.g. [ij, [la, U^j 123, l23j 
|25| . Two underlying ideas are used for these kinds of 
encoder: 

- reducing the coder range M to cover the same area 
around the predicted value for the [k + l)th obser- 
vation y[k -t- 1], y[k -|- 1] G y[k -f 1]. This means 
that the quantizer range M is updated with time 
and a time-varying quantizer (with different values 
of M for each instant, M = M[fc]) is used. Using 
such a "zooming" strategy it is possible to increase 
coder accuracy in the steady-state mode, and, at 
the same time, to prevent coder saturation at the 
beginning of the process; 

- introducing memory into the coder, which makes it 
possible to predict the (A:-|-l)th observation y[A:-t-l] 
with some accuracy and, therefore, to transmit over 
the channel only the encrypted innovation signal. 

In this paper a first-order coder is considered, where 
the predicted value y[k -\- 1] is taken equal to y[k] 
[19I [25I [2a |. In order to describe it, introduce the se- 
quence of central numbers c[k], fc = 0, 1, 2, . . . with initial 
condition c[0] = 0. At step k the coder compares the cur- 
rent measured output y[k] with the number c[A:], forming 
the deviation signal dy[k] — y[k] — c[k]. Then this signal 
is discretized with a given M ~ M[k] according to ([S]). 
The output signal 



dy[k]^q{dy[k],M[k]) 



(6) 



is represented as an i?-bit information symbol from the 
coding alphabet and transmitted over the communication 
channel to the decoder. Then the central number c[fc-|- 1] 
and the range parameter M[k] are renewed based on the 
available information about the master system dynamics. 
Assuming that the master system output y changes at a 
slow rate, i.e. that y\k -\- 1] k, y\k] we use the following 
update algorithms: 



where < p < 1 is the decay parameter. Moo stands 
for the limit value of M[k]. The initial value Mq should 
be large enough to capture the region of possible initial 
values of j/o- Equations (O, ([H]), (jH]) describe the coder al- 
gorithm. A similar algorithm is realized by the decoder. 
Namely, the sequence of M[k] is reproduced at the re- 
ceiver node utilizing ([5]); the values of dy[k] are restored 
with given M[k] from the received codeword; the central 
numbers c[k] are found in the decoder in accordance with 
([7]). Then y[k] is found as a sum c[k] H- dy[k]. 



IV. EVALUATION OF SYNCHRONIZATION 
ERROR 

We now find a relation between the transmission rate 
and the achievable accuracy of the coder-decoder pair, 
assuming that the growth rate of yi{t) is uniformly 
bounded. Obviously, the exact bound Ly for the rate 
of y{t) is Ly — sup|Ca;|, where x is from ^. To 

analyze the coder-decoder accuracy, evaluate the up- 
per bound A = sup|(5y(t)| of the transmission error 

t 
^yi^) — yi(t) ~ yi{t)- Consider the sampling interval 
[tk,tk+i] assuming that dyi[k] < 2M, where dyi[k] = 
yi{tk) — c[k]. Since yi[k] — c[k] + M sign{dyi[k]) it is 
clear that |(5i,(ife)| < M. Additionally, the error may in- 
crease from tk to tk+i due to change of yi(t) by a value 

ik+i 

not exceeding sup \yi{t) - yi{tk)\ < J \yi{T)\dT 
tk<t<tk+i tk 

tk+l 

< J Lydr = LyTs- Therefore the total transmission 

tk 
error for each interval [tk^tk+i] satisfies the inequality: 



\5y{t)\<M- 



LyTs 



(9) 



Inequality ^ shows that in order to meet the inequality 
|(5j,(t)| < A = 2M for all f, the sampling interval Tg 
should satisfy condition 



Ts < A/Ly. 



(10) 



Subtracting Eq. ^ from Eq. ([2]) and taking into 
account the control law ^ we derive an equation for the 
synchronization error in the form 



e{t) =AKeit) +BC{t)-BK5y{t), 



(11) 



c[k + I] ^ c[k] + dy[k], c[0] = 0, A: = 0,1,. 



(7) 



where Ak = A- BKC, C(t) = ip{y2{t))-^{yi{t)). 

Evaluate the total guaranteed synchronization error 
Q — sup lim ||e(i)||, where || • || denotes the Euclidian 

t — *oo 

norm of a vector, and the supremum is taken over all ad- 
missible transmission errors Sy (t) not exceeding the level 
A in absolute value. The ratio Cg = Q/A (the relative 
error) can be interpreted as the norm of the transforma- 
tion from the input function Sy{-) to the output function 



e(-) generated by the system pT|) . Owing to the nonhn- 
earity of the equation (jlip evaluation of the norm Cg is 
nontrivial and it even may be infinite for rapidly growing 
nonlinearities fiy). To obtain a reasonable upper bound 
for Ce we assume that the nonlinearity is Lipschitz con- 
tinuous along all the trajectories of the drive system ([2]). 
More precisely, we assume existence of a positive number 
L^ > such that 

\^{y) - V{y + 5)\ < L^\6\ 

for all y = Cx, x € fl, where il is a set containing all 
the trajectories of the drive system Q, starting from 
the set of initial conditions flo, \S\ < A. For Lipschitz 
nonlinearities ({t) satisfies inequality \C{t)\ < L^\e{t)\. 
After the change K ^ K + L^p, the error equation (fTTj) 
can be represented as 



where 



e{t) = AKe{t) + B^{t) - B{K + L,^)5y{t), 



(12) 



where the variable ^(t) ~ L^e{t) + (, apparently, satisfies 
sector inequality £_{t)e{t) > for all t>0. 

The problem is reduced to quantifying the stability 
properties of ((T^) for bounded input Sy{t). We first an- 
alyze behavior of the system p^ for Sy{t) = 0. To this 
end, we find conditions for the existence of a quadratic 
Lyapunov function V{e) = e^Pe and controller gain K 
satisfying inequality V{e) < —^iV{e) for some /i > for 
dy(t) = and for all ^ satisfying the quadratic inequal- 
ity £,£ > 0. Such conditions can be derived from the 
Passification theorem [23, [2g|, see Appendix. Namely, 
such V and K exist if and only if the transfer func- 
tion of the linear part of the system models H]), ^ 
W{X) — C {XI — A)"^ B is hyper- minimum phase (HMP). 
Recall that the HMP property for a rational function 
W{X) = 6(A)/a(A) where a(A) is a polynomial of de- 
gree n, 6(A) is a polynomial of degree not greater than 
n — 1 means that b{X) is a Hurwi tz p olynomial of degree 
n — 1 with positive coefficients |28j . Now consider the 
case Sy{t) 7^ 0, assuming that the HMP condition holds 
and matrix P and gain K are chosen properly and the 
modified Lyapunov inequality PAk + AJ.P < —fiP is 
valid for some /i > 0. Evaluating the time derivative 
of function V{e) along trajectories of ([2]), ([T]) with ini- 
tial conditions in CIq, using standard quadratic inequality 
\e'^PB\ < y/V{e)y/ViB) after simple algebra we get 

V < -^iV + le^PBiK + L.^)Sy\ < -^iV + Wiy, 



where f = ^/V{B){\K\ + L^)A. Since V < within 
the set ^/V> ii~^v, the value of lim supV{t) cannot ex- 

t — >oo 

ceed A^(L<^ + l-f'^l) Amax(-P)/A*^- In view of positivity of 
P, A„,i„(P)||e(i)||2 < V{t), where A„,i„(P), Xn^^^P) are 
minimum and maximum eigenvalues of P, respectively. 
Hence 



lim|le(i)|l<C+A, 



(13) 



C+ = 



Amin(-P) M 



(14) 



The inequality P^ shows that the total synchroniza- 
tion error is proportional to the upper bound on the 
transmission error A. 

As was shown in the authors' previous paper [l9|, a 
binary coder is optimal in the sense of bit-per-second 
rate, and the optimal sampling time Tg for this coder is 



T, = A/i^Ly), , 



(15) 



where (3 ~ 1.688. Then the channel bit-rate R = l/T^ is 
as follows: 



R = l3Ly/A, 



(16) 



and this bound is tight for the considered class of coders. 
Taking into account the relation P^ for optimal trans- 
mission rate, the synchronization error can be estimated 
as follows: 



lim||e(i)|| <C+pLy/R, 



(17) 



i.e. it can be made arbitrarily small for sufficiently large 
transmission rate R. 

Rem,ark 1. Related estimates for synchronization er- 
rors in coupled systems were obtained in several papers 
m, [33, [m, in, [33. However, in [H [s^, [U, [sl, [s^ either 
existence of Lyapunov functions, i.e. stability of uncou- 
pled systems is required, or a partial stability (stability 
of the synchronization manifold) is provided by a strong 
coupling playing the role of state feedback in the error 
system. In this paper only output feedback is allowed 
and coupling is applied through the control term Bu, i.e. 
in a restrictive manner. That is why the result holds un- 
der the additional assumption (passifiability) caused by 
the nature of controlled problems. Then the partial sta- 
bility conditions are provided by linear observer theory. 
In addition the final result (fT7)) is presented in terms of 
transmission rate, i.e. appeals to the information theory 
view. 

Rem,ark 2. One can pose the following problem: eval- 
uate upper and lower bounds for Ce based on worst case 
inputs Sy{t). Such a problem is similar to the energy con- 
trol problem for systems with dissipation [34, 35] and C'e 
can be interpreted as the excitability index of the system. 
Employing a lower bound for the excitability index for 
passive systems [33, ll^ we conclude that if the gain vec- 
tor K is chosen to ensure strict passivity of the system 
PT|) then the lower bound for Ce is positive, i.e. 



sup lim ||e(t)|| > C^ A. 
|5„(t)l<A*^°° 



(18) 



Thus, for finite channel capacity the guaranteed syn- 
chronization error is not reduced to zero being of the 
same order of magnitude as the transmission error. 



Remark 3. It worth noting that relations ([T3| , p8|) give 
overestimates for the synchronization error, because they 
provide upper and lower bounds for the worst case of the 
input signal Sy (t) . Also the estimate of the mean square 
value of the synchronization error may be used. There is 
a significant body of work in which the quantization error 
signal Sy{t) is modelled as an extra additive white noise. 
This assumption, typical for the digital filterin g th eory, 
is reasonable if the quantizer resolution is high [3^| , but 
it needs modification for the case of a low number of 
quantization levels [T^. The dependence (fT7|) will be 
used for numerical analysis in Sec. [V] 



EXAMPLE. SYNCHRONIZATION OF 
CHAOTIC CHUA SYSTEMS 
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Let us apply the above results to synchronization of 
two chaotic Chua systems coupled via a channel with 
limited capacity. 

Master system. Let the master system ([l} be repre- 
sented as the following Chua system: 



±1 ^p{~xi + ip{yi)+X2), i > 0, 
±2 = a;i — 2:2 + X3 
X3 = -qx2, 
yi{t) ^xi{t), 



(19) 



where yi{t) is the master system output (to be trans- 
mitted over the communication channel), p, q are known 
parameters, x — [xi,X2,X3Y G M"^ is the state vector; 
f{yi) is a piecewise-linear function, having the form: 



ipiy) = moy + mi{\y + l\- \y - 1|), 



(20) 



where niQ, mi are given parameters. 

Slave system. Correspondingly, the slave system equa- 
tions llj for the considered case become 




t >0, 



(21) 



where 2/2(0 i^ the slave system output, z — [zi, Z2, z^]"^ e 
R^ is the state vector, (p{y2) is defined by (PO)) . 
Controller has a form 



u{t) = ~Ke{t), 



(22) 



where e{t) = 2/2 (i) —yiit); yi{t) is a master system out- 
put, restored from the transmitted codeword by the re- 
ceiver at the slave system node (see Fig. [T]); the gain K 
is a design parameter. 

Coding procedure has a form ([S])-®. The input signal 
of the coder is yi (t) . The reference input yi (i) for con- 
troller ([22]) is found by holding the value of yi [k] over the 
sampling interval [kTs, [k + ^)Ts), fc = 0, 1, . . . . 

The following parameter values were taken for the sim- 
ulation: 



FIG. 2: (Color online) Time histories of the state variables of 
master and slave systems (|19|l . (|2ip for A = 1 {Ts = 13 ms, 
R = 75 bit/s): a) xi{t) (dotted line, zi{t) (solid line); b) 2:2 (i) 
(dotted line, 22 (t) (solid line). 
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b) 1, ^w-^w 
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FIG. 3: (Color online) Synchronization error time histories: 
a) e{t) = yi{t) - y2{t), b) e2(i) = X2{t) - Z2{t) for A = 1 
{Ts = 13 ms, i? = 75 bit/s). 
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FIG. 4: (Color online) Time histories of transmission error 
5y{t) (solid line) and coder range parameter M[k] (dashed 
line) for A = 1 (T^ = 13 ms, 7? = 75 bit/s). 
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FIG. 5: (Color online) Time histories of J/i(i) (dashed line), 
yi{t) (solid line) for R = 75 bit/s, t e [48, 50] s. 
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FIG. 6: (Color online) Normalized transmission error Qy vs 
transmission rate R. 
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FIG. 7: (Color online) Normalized state synchronization error 
Q vs transmission rate R. 
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FIG. 8: (Color online) Overall synchronization error gain G 
vs controller gain K. 



Chua system parameters: p = 10, q = 15.6, mo = 
0.33, TOi = 0.945; 

the bound Ly for the rate of yi(t) was evaluated 
by numeric integration of ([T]) over the time interval 
t e [0, ifin], ifin = 1000 s, as Ly = 45; 

parameter A was taken for different simulation runs 
as A = 0.2, 0.4, ...,3.0; 

the sample interval T^ was found for each A from 



- the following coder parameters Mq, A/qo, p in ([5]) 
were taken: AIq — 5, M^o — A/2 (different for each 
A), p = exp(-0.ir,); 

- initial conditions for master and slave systems: 
X, = 0.3, z, = (^ = 1,2,3); 

- the simulation final time tfjn = 1000 s. 

The following accuracy indexes were calculated: 
- relative transmission error 



max 1(5^,(^)1 

0.8tfi„<t<tfi„ ' "' " 



normalized state synchronization error 



(23) 



Q- 



max el 

0.8tti„<t<ttin 



max 

0<t<tti, 



<t)\\ 



(24) 



where Sy{t) = yi{t) — Vi{t), e{t) = x{t) — z{t). Expres- 
sions (j23|) . (|24p characterize transmission error and syn- 
chronization error near a steady-state mode. 



Results of the system examination for K = 1.0 are re- 
flected in Figs. [2HZ1 Figure [H shows time histories of the 
state variables of the master and slave systems (fT9l) , (|2T|) 
Xi{t), zi{t) (plot a) and X2{t), 2:2 (t) (plot b) for A = 1 
{Ts = 13 ms, R = 75 bit/s). The corresponding syn- 
chronization errors ei{t) and 62 (i) are depicted in Fig.[3l 
It is seen that the transient time is about 30 s, which 
conforms with the chosen value of decay parameter p in 
(HI), p = exp(-O.lTs) = 0.9987. Time histories of the 
transmission error Sy and the coder range M[k] are plot- 
ted in Fig. [?] showing that the transmission error actually 
does not exceed M[A:]. Time histories of yi{t) and yi{t) 
for, t G [48, 50] s are depicted in Fig. [5] to provide better 
imagination about the coding procedure. 

Normalized transmission error Qy as a function of the 
transmission rate R is plotted in Fig. [51 Based on the 
simulation results the least-square estimate of the factor 
Gy of the relation between the transmission error Qy and 
the channel rate R having the inversely proportional form 
Qy = Gy/R was numerically found as Gy = 12.2 and the 
curve for Gy/R is also plotted in Fig. [H] (dashed line). 

Synchronization performance may be evaluated based 
on the normalized state synchronization error Q, (|24p . 
shown in Fig. [7] as a function of the transmission rate 
R. The simulation results make it possible to evaluate 
the parameter Gy in the inversely proportional function 
Q = G/R. For the considered example G — 4.0. The 
corresponding curve is plotted in Fig. [7] (dashed line). 
One may notice that the synchronization error is less than 
the transmission error, G < Gy (also compare Figs. [5] 
and [7]). This demonstrates the filtering abilities of the 
synchronization scheme. Based on a theoretical bound, 
simulation data are smoothened with a hyperbolic curve. 
As seen from Fig.|S]it is not necessary to increase control 
power (controller gain) in order to approach minimum of 
limit synchronization error. The value of the gain K does 
not influence synchronization error for sufficiently large 
K. 



VI. CONCLUSIONS 

Limit possibilities of controlled synchronization sys- 
tems under information constraints imposed by limited 
information capacity of the coupling channel are evalu- 
ated. It is shown that the framework proposed in | 19| , is 
suitable not only for observer-based synchronization but 
also for controlled master-slave synchronization via com- 
munication channel with limited information capacity. 

We propose a simple first order coder-decoder scheme 
and provide theoretical analysis for multi-dimensional 
master-slave systems represented in the Luric form. An 
output feedback control law is proposed based on the 
Passification theorem [23, [2g|. It is shown that upper 
and lower bounds on the limit synchronization error are 
proportional to the maximum rate of the coupling signal 
and inversely proportional to the information transmis- 
sion rate (channel capacity). More complex coding pro- 



cedures may provide better synchronization properties. 

The results are applied to controlled synchronization 
of two chaotic Chua systems coupled via a channel with 
limited capacity. Simulation results confirm theoretical 
investigations. 

Future research is aimed at examination of more com- 
plex system configurations, where control signal is sub- 
jected to information constraints too. 
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APPENDIX 



Consider a linear system 

e^ Ae + B(,{t), e = Ce 



(25) 



with transfer function W{X) = G{XI ~ Ay^B= 
6(A)/a(A), where 6(A), a(A) are polynomials, degree of 
a(A) is n, degree of 5(A) is not greater than n — 1. The 
system is called hyper-minimum phase (HMP), if 6(A) is 
Hurwitz polynomial of degree n — \ with positive coef- 
ficients. To find existence conditions for quadratic Lya- 
punov function we need the following result. 

Passification theorem ]2't . [2^]. There exist positive- 
definite matrix P = P'^ > and a number K such that 

PAk + A\P<{), PB = C'^, Ak=A-BKC (26) 

if and only if W{X) is HMP. 

Consider a linear system with feedback 

e = AKe + BS.{t)., e = Ce, Ak = A-BKG. (27) 

Let us show that there exist a quadratic form V{e) = 
e'^ Pe and a number K such that time derivative ^(e) of 
V{e) along trajectories of ([77)1 satisfies relation 



V{e) < for Ce > 0, X 7^ 



(28) 



if and only if W{X) is HMP. To this end assume that K 
is fixed. Relation (|28p is equivalent to existence of the 
matrix P = P"^ > such that e'^ P{AKe + BS)+(,Ce < 
for X ^ 0. Since ^ is arbitrary, the latter in turn, is 
equivalent to matrix relations PAk + A\P < 0, PB = 
G"^ and, by Passification theorem, to HMP condition. 

Remark. It also follows from Passification theorem 
that if HMP condition holds then K satisfying ([^5)1 can 
be chosen sufficiently large. Besides, zero matrix in the 
right hand side of the inequality in (|26p can be replaced 
by matrix — /xP for sufficiently small ^ > 0. 
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